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ABSTRACT 

Planets embedded within debris disks gravitationally perturb nearby dust 
and can create clumpy, azimuthally asymmetric circumstellar ring structures that 
rotate in lock with the planet. The Earth creates one such structure in the solar 
zodiacal dust cloud. In an edge-on system, the dust "clumps" periodically pass in 
front of the star as the planet orbits, occulting and forward-scattering starlight. 
In this paper, we predict the shape and magnitude of the corresponding transit 
signal. To do so, we model the dust distributions of collisional, steady-state 
exozodiacal clouds perturbed by planetary companions. We examine disks with 
dusty ring structures formed by the planet's resonant trapping of in-spiraling 
dust for a range of planet masses and semi-major axes, dust properties, and disk 
masses. We synthesize edge-on images of these models and calculate the transit 
signatures of the resonant ring structures. The transit light curves created by 
dusty resonant ring structures typically exhibit two broad transit minima that 
lead and trail the planetary transit. We find that Jupiter-mass planets embedded 
within disks hundreds of times denser than our zodiacal cloud can create resonant 
ring structures with transit depths up to ~ 1CT 4 , possibly detectable with Kepler. 
Resonant rings produced by planets more or less massive than Jupiter produce 
smaller transit depths. Observations of these transit signals may provide upper 
limits on the degree of asymmetry in exozodiacal clouds. 

Subject headings: circumstellar matter — interplanetary medium — planet-disk 
interactions — methods: numerical — planetary systems 



1. Introduction 



The Kepler Mission is searching for Earth-size planets in edge-on extrasolar plane- 
tary systems by detecting small photometric variations as planets transit in front of their 
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host stars (IBorucki et al.l I2010T ) . To achieve this goal, Kepler was designed to have a 
photometric precision of ~ 2 x 10~ 5 over 6.5 hour timescales for V = 12 Sun-like stars 



Jenkins et al. 


2010; 


Koch et al. 


2010) 



oplanet transits has enabled a wealth of other observatio ns, including aste roseismological 



studies ( iGilliland et al.l 120101). ellipsoida l star variations ( [Welsh et al.l 120101 ). and Doppler 



beaming binaries ( jBloemen et al.l 1201 il l. Kepler m ay also detect photometric va riations 



caused by large satellites orbiti ng transiting planets (|Sartoretti. fc Schneider! Il999l ) as well 



as large, dense planetary rings (IBarnes fc Fortneyll2004j ). 



Recent resolved images of the outer regions of several debris disks hint at another source 
of photometric variability: debris disk structures. These resolved images reveal asymmetric 



from unseen extrasolar plan ets (e.g 



Greaves et al. 


1998; 


Kalas et al. 


2005; 


Golimowski et al 



like ( 3 Pictoris (IHeap et al. 



20061 ; ISchneider et al.l 12009 ) . Asymmetric structures ha ve been observed in ed ge-on disks 



200(1 iGolimowski et all 120061 ) and AU Microscopii dKrist et al. 



20051 ). and othe rs have been observed to slowly rotate on timescales consistent with Keplerian 
orbital motion ( iGreaves et al.l 120051 ) . Several massive hot disks have also been detected in 



the inner regions of nearby stars, where the orbital period is much shorter (e.g. ISong et al. 



20051 ; lAbsil et al.l 120061 ; iBeichman et al.l 120061 ) . Could short-period asymmetries in the inner 
regions of edge-on disks produce photometric variations detectable with Kepler! 

Models of disk-pl anet interactions predi ct that the inner regions of debris disks may 



contain asymmetries (jStark fc Kuchner 



20081 ). As dust grains spiral inward via Poynting- 



Robertson (PR) drag, planets gravitationally perturb their orbits and can trap the dust into 
exterior mean motion resonances (MMRs). This trapping can produce large-scale structure 
in the disk in the form of an overdense, azimuthally asymmetric circumstellar ring that orbits 
in lock with the perturbing planet. 

At least one such structure exists in the inner regions of our own debris disk, the zodiacal 
cloud. Infrared observations of the zodiacal cloud obtained with the DIRBE instrument on- 



board COBE confirmed excess flux leading and trailing the Earth's orbit dReach et al.lll9 95) . 
interpreted as a circumsola r ring of dust resonantly trapped by the Earth (IJackson & Zook 
19891 ; iDermott et al.lll994f ). The inner regions of other debris disks, possibly much more 
dense than the zodiacal cloud, may harbor similar structures. 

A resonant ring structure created by a pl anet on a circular orbit in an exozodiacal cloud 
typically exhibits the following asymmetries (jStark fc Kuchner! 120081 ): 



1. A density deficit, or "gap," near the location of the planet 

2. A density enhancement, or "clump," trailing the planet 
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3. A second clump leading the planet, typically with density less than the trailing clump 

The location of the leading and trailing clumps vary with parameters like planet mass and 
dust size, but for a planet with mass greater th an a few Earth masses, t hese clumps typically 



lead and trail the planet by ~ 90° in its orbit (jStark fc Kuchnerll2008l ) 



For an edge-on system, these density enhancements would produce a variable photo- 
metric signal as the structure rotates in lock with the planet's orbit. As a clump transits, 
the disk blocks additional starlight, producing a stellar extinction light curve with amplitude 
~ Ar r , where Ar r is the change in the radial optical depth along the line of sight to the star. 
This amplitude can be approximated as 

Ar r ~ At A (1) 

where r is the face-on optical depth, R is the radial width of the resonant ring structure, 
and H is the scale height of the resonant ring. For a disk with a range of inclinations on 
the order of 10°, R/2H is on the o rder of a few for a high-c ontrast resonant ring structure 



( jStark fc Kuchnerll2008l ). However, IStark fc Kuchnerl ( 120081 ) showed that At can approach 



lOft, in the absence of grain-grain collisions, where r b is the background optical depth, i.e. 
the optical depth interior or exterior to the resonant ring structure. Exozodiacal clouds 
hundreds to thousands of times more dense than our zodiacal cloud, which has an optical 
depth r ~ 10 -7 near 1 AU, may be capable of producing a signal detectable with Kepler. 

However, another effect works against this signal. The projected angular size of a dust 
clump is typically much larger than the angular size of the star. During the transit of a dust 
clump the majority of the clump, which is near the line of sight to the star but not occulting 
the star, forward scatters starlight toward the observer; forward scattering increases the 
disk flux during dust clump transits. Although this effect is likely small compared to the 
occultation of starlight, a robust calculation of an exozodiacal cloud's light curve must also 
include the disk's scattered light. 

In this paper, we model the transit light curves of collisional exozodiacal clouds with 
asymmetric structures. We limit our investigation to models of ring structures created by 
resonantly trapped in-spiraling dust. We do not model other scenarios that may lead to 
structure, including resonant trapping of planetesimals and transient collisional phenomena, 
but briefly discuss these alternative scenarios in Section 14.11 We show sample light curves 
in Section [3] and determine the disk and planet parameters that would lead to a maximum 
signal amplitude. In Section @] we discuss whether such disks could exist and whether we 
could currently detect their transit light curves. 
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2. Numerical Method 

We modeled the light curves for structured exozodiacal clouds using a two-step process. 
First, we modeled the 3D dust distributions for steady-state collisional debris disks perturbed 
by single planet orbiting a Sun-like star. Second, we illuminated these models with starlight 
and synthesized edge-on images as a function of planetary phase. 



2.1. Exozodi Models 

For a signal ~ 10~ 4 , possibly detectable by Kepler, the disk's background optical depth 
lb ~ 1CT 5 , i.e. hundreds of times greater than the zodiacal cloud. In such dense disks the 
average time for two dust grains to collide, £ co n, can be shorter than the Poynting- Robert son 
time, tpR, for the dust grains that dominate the disk's optical depth. Modeling the light curve 
for these disks therefore requires a model that treats both resonant gravitational dynamics 
as well as grain-grain collisions. 

We produced models of collisional exozodiacal clouds with resonant ring structures us- 
ing a collisional grooming algorithm. This algorithm simultaneously solves the equations of 
motion governing small dust grains in a planetary system and the number flux equation in- 
cluding destruction of dust grains via grain-grain collisions. In short, the collisional grooming 
algorithm takes a collisionless debris disk "seed model" as input and iteratively "grooms" the 
disk until the 3D dust distribution matches that of a steady-state collisional disk with a given 
dust produc tion rate. A more detaile d explanation of the collisional grooming algorithm can 



be found in IStark &: Kuchnerl ( 120091 ) . 



To produce the collisionless seed models, we simulated the dynamical interactions be- 
tween a cloud of dust and a single planet on a circular orbit around a Sun-like star. We 
investigated four different planet masses, M p = [1 M E , 5 M E , 1 M N , 1 Mj], where M E , M N , 
and Mj are the mass of the Earth, Neptune, and Jupiter, respectively. Each of these planet 
masses were modeled with two values for orbital semi-major axis, a p = [0.5, 1.0] AU. For the 
Jupiter-mass case we also created models with a p = 0.1 AU. 

For each of these models, we integrated the orbits of 12,500 dust grains launched from 
a belt of parent bodies exterior to the orbit of the planet. We used a hybrid symplectic 
i ntegrator with an integ ration step size equal to one hundredth of the period of the planet 



(IStark fc Kuchnerl 120081 ). The parent body orbits were distributed uniformly in semi-major 
axis between 2.5a p and 3.0a p , uniformly in eccentricity between and 0.1, and uniformly in 
inclination between and 10°. All other parent body orbital parameters were distributed 
uniformly between and 27r. 
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For the seed model integration, we distributed the dust grains in equal numbers among 
25 logarithmically-spaced values, ranging from 0.43355 to 0.00046, where /3 is the ratio 
of the radiation pressure force on a grain to the stellar gravitational force. Each discrete 
value represents the behavior of a small range of (3 values from /3/\/L33 to (3 a/1-33, so 
that we represent grains down to the blowout size (/3 = 0.5). These (3 values correspond 
to 25 values of grain size, ranging roughly from sub-micron to nearly millimeter sizes. We 
calculated s((3) according to 

n (Qpr) 1 g cm~ 3 

s = 0.57 — jum, (2) 

P P 

valid for a Sun-like star, where p is the mass density of a dust grai n and (Qpr) is th e 



Poynting- Robert son efficiency, Qpr, averaged over the solar spectrum (jBurns et al.lll979l ). 
We let p — 2 g cm -3 , appropriate for "rocky" grains. Qpr, an implicit function of s, is given 
by 

QPR = Qabs + Qsca (1 - (cos a)) , (3) 

where Q a b s and Q sca are the absorption and scattering efficiencies of the dust grain and 
(cos a) is the cosine of the scattering angle w eighted by the scattering phase function and 



averaged over all directions (jBurns et 



an gle w e] 
~all ll979h . 



Values for Q abs and Q sca , as well as the scattering phase function for exozodiacal dust 
grains are unknown. Observations of dust in the outer regio ns of debris disks reveal grains 



that appear to deviate significantly from Mie theory (e.g. iKalas et al.l 120051 ; iKrist et al. 



20101 ). In light of this, we chose simple generic laws for these quantities and investigated the 
effects of varying the dust grains' scattering parameters. We let Q a b s = 1 for A < 2ns and 
Qabs = (2ns/\) 2 for A > 2ns. We used three values for the albedo, u = [0.1,0.3,0.5], where 
y = Q sna /(Q sr , a + Qah s). Finally, we used a Henyey-Greenstein scattering phase function 



(IHenyey fc Greensteinlll94ll ). given by 



P (") = — 777 , (4) 

47r(l + ^_2^cosa) 3/2 

and investigated three values for the scattering asymmetry parameter g = [0.2,0.5,0.8]. 
Since g = (cos a) , g — corresponds to isotropic scattering and increasing values of g in the 
range < g < 1 correspond to increasing degrees of forward scattering. 

We recorded the positions and velocities of the grains in the frame co-rotating with the 
planet every £pr/14300 years for a maximum time of 1 Gyr. This time between records 
resolves the collision time for the smallest grains in our 500 zodi simulations by a factor 
of ~ 10 4 . During the integration, grains were removed if their semi-major axis a < a p /3 
or a > 300 AU, if they collided with the planet given realistic planet radii, or if they were 
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scattered into hyperbolic orbits and ejected from the system. For the case of M p = 5 
Me, we set t he planet radius eq ual to 1.5 Earth radii, appropriate for a roughly Earth-like 



composition (ISeager et al.l 120071 ). For all other planets we used the radii of their respective 



solar system counterparts. 

We then processed the collisionless seed models with the collisional grooming algorithm 
to obtain steady-state collisional disk models. We weighted the dust production rate of 
each grain size such that at the moment of launch from a planetesi mal the dust grain s were 



produced according to a Dohnanyi crushing law, dN/ds oc s -3 5 ( iDohnanyil Il969t ) . The 
collisional grooming algorithm then solved the number flux equation and self-consistently 
determined the size distribution at all points in the disk. 

For each collisionless seed model, we produced 6 different collisional disk models with 
increasing optical depths (i.e. increasing disk masses/dust production rates). We chose the 
total dust production rates such that the 6 collisional disk models had maximum optical 
depths of 10~ 7 , 5 x 10~ 7 , 10~ 6 , 2 x 10 -6 , 5 x 10~ 6 , and 10" 5 . For the Jupiter-mass case 
we also modeled a maximum optical depth of 5 x 10~ 5 . These maximum optical depths 
roughly correspond to dust production rates of 2.6 x 10~ 10 , 2.3 x 10~ 9 , 7.1 x 10~ 9 , 2.1 x 10~ 8 , 
8.5 x 10 -8 , 2.2 x 10~ 7 , and 2.0 x 10~ 6 Me Myr -1 for dust grain sizes ranging from the blowout 
size up to ~1 mm in radius, though dynamical differences can cause the dust production rate 
to vary by as much as a factor of 2 among models with the same optical depth. Here "optical 
depth" refers to the dust grain cross section per unit area of the disk integrated along a line 
normal to the disk plane and does not directly depend on the dust grains' optical constants. 
To mitigate the effects of Poisson noise when determining the maximum optical depth, we 
azimuthally averaged the disks' optical depths into radial bins of width 0.05a p . For the rest 
of this paper, we will use 1 "zodi" to refer to a maximum optical depth of 10~ 7 ; our models 
range from 1 zodi to 500 zodis of dust. 

We made sure that the collisional grooming algorithm resolved any fine structure within 
the disk by adopting a 512 x 512 x 128 grid with a bin size of 0.05a p x 0.05a p x 0.025a p . 
We ran the collisional grooming algorithm until no single dus t grain record changed by more 



than 5% from one iteration to the next. I St ark Kuchnerl (120091 ) performed a number of 



tests showing that this criterion is sufficient to achieve convergence to the solution of the 
number flux equation for disks with and without azimuthal asymmetries. This degree of 
convergence is roughly equivalent to an uncertainty of less than 0.05% in the bin-averaged 
density distribution for the densest regions of our models. 



Similar to the collisional models of the Kuiper Belt dust disk in lKuchner fc Stark! ( 120101 ). 
our models draw a distinction between cratering events and catastrophic collisions. Catas- 
trophic collisions occur when the collisional energy measured in the center of momentum 
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frame of the two dust grains exceeds the critical disruption energy of the dust grains. In this 
case, both grains are destroyed. All other events may be deemed cratering events, which we 
ignore. In these cases, both grains remain unchanged. In the strength regime the critical 
disruption energy is commonly approximated by 



Q = A S 



1 m 



We used the parameters for "rocky" grains from lKrivov et al.l (120061 ) and let b s 
A s = 10 6 erg g^ 1 . 



(5) 



-0.24 and 



We note that our treatment of collisions ignores fragments. Although this may be a 
decent appr oximation for icy grains produced in the outer regions of the solar system as 
discussed in iKuchner fc Stark! (120 lOf ). the validity of this approximation is not as clear for 
the grains simulated here. We discuss this approximation further in Section 14.31 



2.2. Synthesizing Exozodi Light Curves 

For each of the above models, we used our publicly-available IDL software package 
dustmap to synthesize edge-on images of the dust disk as a function of planetary phase. Given 
a user-specified disk orientation and distance, dustmap illuminates the disk with starlight 
and calculates the scattered and emitted light from each individual dust grain in the disk 
model. The code is capable of calculating the scattering and absorption efficiencies of the 
dust as a function of composition, grain size, and wavelength using Mie theory, but we chose 
to use the generic values described in Section 12.11 For thermal emission images, dustmap 
uses the stellar spectrum, approximated as a blackbody, and the wavelength-dependent dust 
absorption efficiency to calculate the equilibrium temperature of the dust as a function of 
grain size and circumstellar distance. For scattered light images, dustmap determines the 
scattering angle for each dust grain and uses a Henyey-Greenstein phase function to calculate 
the flux scattered toward the observer. 

When calculating the contribution of the stellar flux, our imaging code takes into ac- 
count extinction of the stellar light by transiting dust, but does not include the effects of 
stellar limb darkening. For transits of optically-thick objects with well-defined edges, such 
as planets, these effects are important. However, the clumps in an optically thin exozodiacal 
cloud have comparatively smooth features; the projected angular size of the dust clump's 
"edge" is typically larger than the angular size of the star, rendering limb darkening effects 
unimportant. 

We imaged our model disks at a wavelength of 575 nm, the peak in Kepler's response 
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curve (IKoch et al.ll2010l ). We use our images at this single wavelength to approximate the 
exozodi signal integrated across the entire Kepler photometric response curve, which extends 
roughly from 425 nm to 900 nm. This approximation is valid given our approximations for 
Qabs and Q 

sea' 

We produced transit light curves from our models by calculating the total flux at 72 
values of planetary phase equally spaced over 360°. The relative intensity of our transit light 
curves was calculated as 

"max 

where F(4>) is the total (disk plus stellar) flux as a function of planetary phase 0, F max is the 
maximum of F(<f>), and the stellar flux was calculated assuming a 1 solar radius spherical 
blackbody with an effective temperature of 5778 K. 



3. Results 

3.1. Resonant rings in collisional disks 

Figure [TJ shows face-on optical depth histograms for twelve of our models. Each row 
corresponds to one of the four modeled planet masses and each column corresponds to one 
of the seven modeled zodi levels. Each model is scaled independently to show detail. The 
location of the star and planet are marked with a star and a circle, respectively, and the 
planets orbit counterclockwise. The birth ring, which traces the location of the planetesimals, 
appears as a symmetric circular ring between 2.5 and 3.0 AU. Resonant ring structures appear 
interior to the birth ring, near 1 AU. 

The 1 zodi models (left column) illustrate the features typically attributed to resonant 
ring structures: a gap near the location of the planet, an overdense "clump" trailing the 
planet, and in some cases a slightly less overdense clump leading the planet in its orbit. The 
resonant ring created by a Jupiter-mass planet is clearly more dense and asymmetric than 
the resonant ring created by an Earth-mass planet. This is because the Jupiter- mass planet 
traps dust more efficiently, thereby trapping more dust and trapping a larger fraction of that 
dust in the 2 : 1 resonance, which contributes strongly to the structure's bilobed appearance. 

The rows in Figure J illustrate the ef f ects o f colli sions on resonant r i ng; str uctures. As 
previously described in IStark &: Kuchnerl (120091 ) and iKuchner &: Stark! (120101 ). increasing 
the optical depth of a disk increases the collision rate, which in turn erases resonant ring 
structure. However, the top-right histogram in Figured] shows that even for an optical depth 
of 100 zodis, a Jupiter-mass planet can create significant resonant asymmetry in a disk. 
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1 zodi 10 zodis 100 zodis 




-2 2 -2 2 -2 2 

x(AU) 



Fig. 1. — Face-on optical depths for twelve of our models. Each optical depth histogram is 
scaled independently to show detail. For all models shown, u = 0.3, g = 0.5, and a p = 1.0 
AU. Each row corresponds to one of the four modeled planet masses, as labeled along the 
right side. Each column corresponds to one of the maximum optical depths modeled, as 
labeled along the top of the figure. At optical depths ~ 100 zodis, a Jupiter-mass planet can 
create significant resonant asymmetry while an Earth- mass planet cannot. 
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3.2. Light curves 



Figure [2] shows a sample light curve for a 10 zodi disk perturbed by a Jupiter-mass 
planet. The planetary transit occurs at a planet phase of 0°, but is not shown in this plot 
because it is 3 orders of magnitude larger than the disk transit and lasts for only a few 
hours. The trailing and leading dust clumps in the disk transit at approximately 90° and 
270°, respectively. The solid line in Figure [2] shows the total light curve for the star/disk 
system, the result of a competition between stellar extinction and forward scattering from 
the disk. 

The dashed line shows the variation due to stellar extinction alone. Similar to a plane- 
tary transit, the enhanced column density of dust between the observer and the star during 
the transit of a dust clump reduces the amount of starlight that reaches the observer. How- 
ever, unlike a planetary transit, the projected angular size of the light-blocking material 
(the clump) and the angular scale of the clump's "edges" are typically much larger than the 
angular size of the star. This difference in angular size produces transit minima that are 
much broader than that of a planet. 

The dotted line in Figure |2] show the variations in the disk flux alone. Since the angular 
size of a clump is much larger than the angular size of the star, dust just pre- or post-transit 
forward-scatters additional starlight toward the observer. This effect produces a disk flux 
that peaks when the stellar extinction reaches a minimum. 

The forward scattering effect for our g = 0.8 model, shown by the dotted line in Figure 
[21 may be unrealistically large. The outer regions of observed debris disks exhibit q values 



typically much smaller, on the orde r of 0.2 o r less ( lAugereau et al.l Il999l ; iKalas et al.l 12005 



Debes et al.l 120081 ; iKrist et al.l 120101 ) . iHongl (119851 ) fit the scattering phase function of the 



zodiacal cloud using a weighted series of three Henyey-Greenstein phase functions and found 
the dominant term to be highly forward scattering, with g = 0.7. However, this assumes 
a dust distri b ution oc r -1 , and when correcting for the r -L34 distribution determined by 



Kelsall et al.l ( 119981 ). the effective g value decreases significantly. 



We produced other models with lesser degrees of forward scattering (g = 0.2, 0.5). 
Forward-scattered light from the disk contributes negligibly to the light curves for all of 
these models. Figure |3] shows an example light curve for such a model. This model is the 
same as that shown in Figure [U except g — 0.2 and the disk density was increased to 100 
zodis. 
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Fig. 2. — Light curve for a 10 zodi dust disk with oj = 0.5 and g = 0.8, perturbed by a 
Jupiter-mass planet on a circular orbit with a p = 0.1 AU. Transiting dust clumps create the 
broad minima. The total light curve is the result of a competition between occulted stellar 
flux and forward-scattered flux from the disk. 
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Fig. 3. — Light curve for a 100 zodi disk with g = 0.2 and all other parameters the same as 
in Figure [2j The transit depth amplitude is approximately 7x 10~ 5 . The trailing and leading 
clumps have similar optical depths, a result of small grains being preferentially trapped to 
librate about the leading center of the 2 : 1 resonance. 
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3.3. Asymmetry reversal and light curve complexities 

The light curves shown in Figure [2] and Figure [3] also highlight the morphological effects 
of collisions on resonant ring structures. The 10 zodi model shown in Figure [2] shows a 
trailing clump that is significantly more dense than the leading clump, while the 100 zodi 
model shown in Figure [3] shows a trailing clump density on par with the leading clump 
density. At a p = 1.0 AU, the changes are even more significant, as shown in the upper row 
of Figure [TJ In this case, the leading clump is slightly more dense than the trailing clump 
for a disk density equal to 100 zodis. 

It appears that collisions are "choosing" to remove grains from the trailing clump, 
but a more detailed look at the optical depth as a function of grain size reveals that this 
effect stems from resonant dust grain size-sorting. Figure H] shows the optical depth for all 
dust grains (resonantly trapped and non-trapped) within the region of the 2 : 1 MMR for 
the case of a Jupiter-mass planet orbiting at 1.0 AU in a 100 zodi dust disk. These dust 
grains were selected by their semi-major axes a such that 0.95 a 2; i < a < 1.05 a 2:1 where 
0-2:1 = a p 2 2//3 (l — /3) 1 / 3 . The largest grains (/3 < 0.14) create a dust clump that trails the 
planet while the smallest grains (f3 > 0.14) create a dust clump that leads the planet. Figure 
H] shows that for grains within the region of the 2 : 1 MMR, the leading dust clump contributes 
roughly 50% more to the optical depth than the trailing dust clump. Including grains at all 
semi-major axes reduces this asymmetry reversal to approximately 5%. A close examination 
of the paths of individual particles reveals that this asymmetry reversal does not suffer from 
small number statistics; roughly 350 dust grains with > 0.14 occupy orbits that contribute 
to a denser leading clump while 270 contribute to a denser trailing clump. 

This size-dependent asymmetry extends to disks with much lower optical depths, ~ 1 
zodi. To first order, collisions simply select which grain size dominates the optical depth of 
the resonant ring structure. For lower disk optical depths ~ 1 zodi, the PR times of grains 
with (3 < 0.14 can be less than their collision times. In this case, large grains migrate inward 
and reach the planet's MMRs to create a trailing clump that dominates the resonant ring's 
optical depth. For denser disks, like that shown in Figure HJ the grain-grain collision time is 
reduced and only grains with (3 > 0.14 can migrate inward to populate the planet's MMRs, 
creating a dominant leading clump. 

The leading-trailing clump asymmetry only reverses for the M p = 1 Mj models. These 
models are the only models for which j : 1 exterior resonances dominate the MMR distri- 
bution, and the j : 1 resonances are unique in that they can lead to asymmetric libration. 
Asymmetric libration arises from a balance between the direct and indirect terms of the 
disturbing function and splits the libration center of a resonance into two separate libration 
centers for particles that obtain moderate orbital eccentricities. These libration centers lead 
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Fig. 4. — Face-on optical depth of all grains with semi-major axes near the semi-major axis 
region of the 2 : 1 MMR for a 100 zodi simulation with M p = 1 Mj and a p = 1.0 AU. The 
planet, marked with a circle, orbits counterclockwise. The smallest grains (f3 > 0.14) create 
a clump that leads the planet in its orbit while the largest grains (f3 < 0.14) create a clump 
that trails the planet. In this case, the small grains dominate the optical depth of the disk, 
leading to asymmetry reversal with a denser leading clump. 



- 15 - 



and trail the planet by ~ 90° for the M p = 1 Mj case. It is useful to refer to a single j : 1 
MMR with asymmetric libration as two separate resonances, which we will refer to as leading 
and trailing resonances. 

The leading and trailing j : 1 resonances have different trapping efficiencies. Several 
authors have found that in the case of a planet migrating outward toward a collection 
of massless particles, particles trapped in the 2 : 1 resonance typically populate the lead- 
ing and trailing resonances equally or preferentially p opulate the trailing resonance, creat- 
ing a dominant clump of particles trailing the planet (jWyattl 120031 ; iMurray-Clay &: Chiang 
20051 ). However, in their simulations of the outward migration of Neptune in the Solar 
System, IMurray-Clay fc Chiang! (120051 ) showed that for some migration rates the particles 
preferentially populated the leading resonance by up to a factor of 2 (see their Figures 
15 and 16), creating a n asym metry reversal similar to what we observe in our models. 
Murray-Clay fc Chiang! ( 120051 ) attributed this asymmetry reversal to an asymmetric shift 
in the stable and unstable equilibrium points for the resonant angle such that they form a 
relatively flat potential near the leading center. As a result, the particles spend more time 
traversing the region near the leading center where the librat ion rate is much lower, leadin g 
to a higher probability of capture into the leading resonance (IMurray-Clay fc Chiang! 120051 ). 
Although our models feature particles spiraling in toward a non-migrating plan et, the physics 



of asymmetric libration is similar to the case of outward planetary migration (lQuillenll2006 



by 



Mustill &: Wvattll201l!) and the o rigin of the asymmetry is likely identical to that described 



Murrav-Clav fc Chiang! (l2005h . 



The asymmetry reversal seen in the j : 1 resonances demonstrates the size-sorting effects 
of resonances. For example, for the case of a Jupiter-mass planet, we expect to find smaller 
grains in the leading clump and larger grains in the trailing clump. These clumps may 
therefore exhibit differing spectral features, and differences in their color may possibly be 
detectable with future multi-band resolved imaging or multi-band transit photometry. 

The transit light curves of resonantly trapped dust are not always composed of two 
distinct minima. Figure [5] shows the light curve for a disk with identical parameters to 
that shown in Figure [31 but perturbed by a 5.0 Me planet on a 0.5 AU orbit. The ~ 10~ 7 
variations in the light curve are attributable to Poisson noise from the model, but larger 
variations are real. In this case, the resonant ring structure does not resemble a simple 
bilobed structure (see Figured]); the multiple, subtle clumps in the resonant ring produce 
several smaller minima. 
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Fig. 5. — Light curve for a 10 zodi disk with u = 0.5 and g = 0.2 perturbed by a 5.0 M E 
planet on a circular orbit at 0.5 AU. The minimum near 180° is real, created by a density 
enhancement in the resonant ring. Variations of amplitude ~ 10~ 7 can be attributed to 
Poisson noise. 
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3.4. Light curve amplitudes 

We would like to know the maximum signal amplitude that exozodiacal resonant ring 
structures can create. To determine this, we calculated the amplitude A of each modeled 
light curve, given simply by A = (1.0 — I m i n )- 

As alluded to by Figure [5j some modeled light curves suffer significantly from Poisson 
noise. To avoid amplitudes that were dominated by Poisson noise, we fit each light curve 
with an eighth-order polynomial and then calculated the residual standard deviation. We 
rejected curves with amplitudes less than four times the residual standard deviation. We 
re-synthesized the rejected curves using a finer sampling for the planet phase (2000 values 
for planet phase as opposed to the original 72) and then binned the light curves down to 
50 planet phase values to mitigate Poisson noise. We fit these curves again with an eighth 
order polynomial, and permanently rejected those with amplitudes less than four times the 
residual standard deviation. In all cases, the final rejected light curves failed our Poisson 
noise test because relatively little resonant ring structure existed within the disk. 

Figure [6] shows the results for the a p = 0.5 AU case with u = 0.5 and g = 0.2, although 
the results are nearly independent of g. Only two data points are shown for the M p = 1 

case. The amplitudes at larger disk densities did not meet our Poisson noise criteria 
described above; an Earth-mass planet creates a very weak resonant ring structure in disks 
with optical depths greater than a few tens of zodis. 

As expected, the light curve amplitude increases with disk density. However, the am- 
plitude of the signal does not scale linearly with the number of zodis; for the 1 Mj curve, 
the maximum signal amplitude for the 100 zodi case is only ~40 times greater than the 1 
zodi case. We ascribe this non-linearity to two effects brought about by grain-grain colli- 
sions. First, the enhanced collision rate in denser disks will reduce the fraction of dust that 



migrates inward from the planet esimal belt to populate the planet's MMRs (jWyattl 12006 



Stark fc Kuchnerl 120091 ) . Second. IStark fc Kuchnerl (120091 ) showed that grain-grain collisions 



preferentially remove dust from MMRs and smooth out disk asymmetries. These phenomena 
effectively curb the signal amplitude as the disk density increases. 

The only exception to this rule is the light curve amplitude for a resonant ring created 
by an Earth- mass planet, which decreases with increasing optical depth. In general, trapping 



efficiency decreases as the dust grain size decreases (e.g. IWyattl 120031 ) . and increas ing the 



optical depth reduces the typ ical dust grain size that reaches the planet's resonances ( jWyatt 



20061 ; IStark fc Kuchnerl 120091 ). In the Earth-mass case, the planet traps dust so inefficiently 
that even a moderate optical depth results in grains that are typically too small for the 
planet to trap. The bottom row in Figure [TJ shows clearly that an Earth-mass planet traps 
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Fig. 6. — Light curve amplitude as a function of zodi level for all four modeled planet masses. 
In all cases a p = 0.5 AU, u = 0.5, and g = 0.2. Jupiter-mass planets can create resonant 
disk structures with maximum transit depths ~ 10~ 4 while Earth-mass planets cannot create 
detectable resonant disk structures. 
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relatively little dust into resonance for optical depths greater than a few zodis. 

Figure M shows that disks perturbed by planets less massive than Neptune are unlikely 
to produce a signal detectable with Kepler. Disks perturbed by Neptune-mass planets can 
produce signals greater than 10~ 5 in amplitude, while dust trapped by super-Earths can 
only exceed the 10 -6 level. Resonant rings created by Earth-mass planets have little hope of 
being detected via transit photometry, with signal amplitudes on the order of 10 -7 or less. 

Figure E] also shows that exozodiacal dust resonantly trapped by a single Jupiter-mass 
planet on a circular orbit can produce a transit signal that exceeds the 10~ 4 level for disks 
with optical depths of several hundred zodis. To within a factor of a few, this was the largest 
amplitude we could produce with our models. We found that even larger zodi values could 
increase the scattered disk flux signal by a factor of a few, but disks with optical depths 
greater than ~ 1000 zodis at a few tenths of an AU have shorter lifetimes and are therefore 
less likely to be detected (see discussion of disk lifetimes in Section BZQ) . 

We briefly examined the effects of increasing the planet mass beyond 1 Jupiter mass 
to M p = 2 Mj and 5 Mj, and found that the signal amplitude decreased. By examining 
the semi-major axis distribution of dust grains, we determined two causes for this decrease. 
First, increasing the planet mass beyond a Jupiter mass improves the trapping efficiency of 
higher-order MMRs (e.g. 3 : 1 and 5:2). Since resonant trapping pumps up the dust grains' 
orbital eccentricities, the increased population of grains in higher-order MMRs serves to 
enhance the collision rate at larger circumstellar distances. This in turn reduces the number 
of grains available to migrate inward to the 2 : 1 resonance, where trapping is most efficient. 

Second, increasing the planet mass moves the exterior chaotic zone boundary to larger 
semi-major axes. The exterior chaotic zone boundary, expressed as 



Ochaotic =a p [l + K ^ ~ - - j J , (7) 

is a semi-major axis boundary int erior to which MMRs overlap and sustained resonant 



trapping is unlikely ( IWisdomlll980f ). Estimates for the value of k vary; we chose the most 



recent estimate, k = 2-0, w hich comes from numerical simulations of the Fomalhaut debris 



ring by IChiang et al.l (120091 ) . If the semi-major axis of an exterior j : k MMR, given by 



/ ■ \ 2/3 

a r . k = a p r-\ (1-/3) 1/3 , (8) 

is less than a c haotic, then the resonance will not be significantly populated. Figure [7] shows 
the semi-major axes for the 2 : 1 and 3:2 MMRs as functions of /3, along with the chaotic 
zone boundaries for several planet masses. As Figure [7] shows, for M p = 1 Mj the 2 : 1 MMR 
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is exterior to the chaotic zone boundary for f3 < 0.4, but for M p = 5 Mj only grains with 
(3 < 0.2 can be trapped in the 2 : 1 MMR, significantly reducing the population of the 2 : 1 
MMR in the latter simulation. 

We also examined the amplitudes of the variable disk flux signal, shown as dotted lines 
in Figures |2j [3j and |5j Disks inclined from an edge-on orientation by an amount greater 
than their disk opening angle would not exhibit the broad transit minima shown in these 
figures, but would still feature the variable disk flux caused by forward-scattering. Our 
models predict that the largest possible forward-scattered light signal that a disk with a ring 
of resonantly-trapped dust can produce is roughly 5 x 10~ 6 . Detecting resonantly-trapped 
dust rings in disks that are not in an edge-on orientation would therefore require photometric 
precision better than one part in 10 6 over the course of an orbital period. Additionally, to 
produce such a large scattered light signal the disk must have an optical depth of at least 100 
zodis, strongly forward scattering dust (g = 0.8) with a large albedo (u = 0.5), a resonant 
ring structure created by a Jupiter- mass planet just exterior to the sublimation distance for 
silicate grains, and an inclination just greater than the disk opening angle. Simultaneously 
satisfying all of these criteria may only occur very rarely and we consider this scenario 
unlikely. 



4. Discussion 

4.1. Disk lifetimes 

Figure [6] shows that the amplitude of a resonant dust ring's light curve approaches the 
10~ 4 level only for massive perturbing planets, ~ 1 Mj, with disk optical depths ~ 100 zodis. 
Additionally, the likelihood of photometrically detecting transiting dust clumps improves as 
the planet semi-major axis decreases, since the transit probability increases, the number of 
orbit-foldings of the light curve increases for a given total observation time, and the required 
timescale for photometric stability decreases. Radial velocity and transit measurements 
have already shown a wealth of close-in massive planets, but should we also expect to find 
hundreds of zodis of dust? 



Mass ive amounts of hot dust (> 1000 zodis 



have already been observed around other 



(e.g. 


Beichman et a 


.2005 


Song et al 


2005 


Wya 


;t et al. 


2005: 


Absil et al.ll2006: 


Rhee et al. 


Akeson et al. 


2009; 


Absil et al.ll2009: 


Moor et al. 


2009; 


Melis et al. 


2010; 


Weinberger et al. 


Millan-Gabet et al. 


2011 


, although the frequency of these disks around mature (> 1 



2008 



2011 



Gyr) solar- type stars is only on the order of a few percent (jBeichman et al.ll2006l ). Obser- 
vations have also revealed a few less massive exozodiacal disks with optical depths greater 




Fig. 7. — Chaotic zone boundaries (solid lines) as functions of (3 for M p — 1, 2, and 5 Mj, 
plotted with the locations of the 2 : 1 and 3:2 MMRs (dotted lines). Sustained resonant trap- 
ping can occur when the MMR semi-major axis is greater than the chaotic zone boundary. 
As planet mass increases the 2 : 1 resonance intersects the chaotic zone boundary at smaller 
values of j3; more massive planets trap fewer small grains into the 2 : 1 resonance. 
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than a few hundred zodis (IMoor et al.ll2009l ; iMillan-Gabet et al.ll201ll ). The least massive of 
these detected exozo diacal clouds are consistent with the long term collisional erosion of a 
belt of planetesimals (IWyatt et al.ll2007l ; IMoor et al.ll2009f ). similar to the disks we modeled. 
But the dust in many of these systems is likely transient in nature; current estimates restrict 
the typical disk mass e- folding times to less than ~ 10 Myr, much less than the age of many 
of these systems (IWyatt et al.l 120071 ). 



However, transient collisions may also produce transient resonant rings. After all, the 
fundamental mechanism behind our models is very basic: collisions produce dust, some 
fraction of which migrates inward and becomes trapped in a planet's mean motion resonances. 
We suspect that a massive collisional avalanche initiated by a single event ends with the 
inward migration of any remaining bound grains. Future modeling of massive, transient 
collisions is necessary to understand the likelihood of this scenario. 



The lifetime of dense exozodiacal disks is not well known, bu tlWvatt et a 



loose upper limits on the age of these disks. Using Equation 21 in lWyatt et al 



. (2007) place 



(120071 ). along 



wi th the fact that the disk lifetime could be two orders of magnitude longer (see Section 3.2 



m 



Wvatt et al.l (120071 )). we calculate maximum disk lifetimes of 1-100 Myr, 34-3400 Myr, 
and 0.17-17 Gyr for our a p = 0.1, 0.5, and 1.0 AU 100 zodi disk models, respectively. These 
calculated lifetimes vary with a p and suggest that disks with 100 zodis of dust at distances 
less than ~1 AU only exist around the youngest stars, while disks with 100 zodis of dust at 
distances larger than ~1 AU may exist for up to billions of years. However, the lifetime varies 
with a p only because our models assumed planetesimals distributed from 2.5a p to 3.0a p ; for 
the case of a p = 0.1 AU the assumed planetesimals are on short-period orbits near 0.3 AU. 
In reality, dust may drift inward to 0.1 AU from a belt of planetesimals located farther out 
where the disk is significantly longer-lived. 

The model shown in the lower right panel of Figure [T] clearly illustrates how a more 
distant source of planetesimals may deliver copious amounts of dust to the very inner regions 
of a planetary system. In this model planetesimals near 3 AU produce 100 zodis of dust. The 
dust migrates inward past the Earth-mass planet at 1 AU, which has little effect on the dust 
distribution. The optical depth interior to 1 AU (but exterior to our integration-imposed 
inner cutoff of 0.3 AU) is approximately constant as a function of circumstellar distance; 
nearly all of the dust that migrates to 1 AU would continue to migrate to 0.1 AU. 

Dense exozodiacal clouds may be fed by even more distant sources of dust, e.g. mas- 
sive Kuiper Belt analog s, fo r which the d i sk's l ifetime would be of no concern. Both 
Kuchner fc Stark! ( 120101 ) and IVitense et al.l (120101 ) simulated the distrib ution of dust pro- 
duced by t he Kuiper Belt using mo dels that included grain-grain collisions. Kuchner &: Stark 
( 120101 ) and lVitense et al.l (120 lOf ) both showed that small Kuiper Belt dust grains can migrate 
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into the inner solar system via PR drag, eve n if the Kuiper B e lt dus t disk is dense enough 
to be in the "collision-dominated" regime. iKuchner fc St arid ( 120101 ) estimated that ~ 60 
zodis of dust would m ake it to 1 AU from a 10 00- zodi Kuiper Belt analog in the absence 
of Saturn and Jupiter. IStrubbe &: Chiang! (120061 ) showed that strong stellar winds can also 
transport large amounts of dust into the inner disk in spite of grain-grain collisions; an AU 
Mic disk-analog could deliver more than a thousand zodis of dust to the inner few AU if the 
stellar mass loss rate were greater than 100 times solar. 

Dust from a more distant source must migrate over a larger distance to reach the inner 
disk and this migration must be completed in an amount of time roughly equal to the grain- 
grain collision time. This means that for a given collision time in a disk, the dominant dust 
grains that arrive from a more distant source must migrate more quickly, i.e. all other things 
being equal, more distant sources will deliver smaller dust grains to the inner disk. Since 
planets with M p < Mn don't trap the smallest dust grains with 100% efficiency, we expect 
the amplitude of the dust transit signal to decrease for these models. More massive planets 
still trap small grains with very high efficiency, though, so we expect our M p = 1 Mj to 
remain relatively unchanged in amplitude for more distant sources delivering equal amounts 
of dust to the inner disk. However, the shape of the transit signal may change slightly, since 
small grains exhibit a stronger asymmetry reversal in our models (see Section [3731) . 



Our models did not consider scenarios in which the planetesimals themselves were 
trapped in external MMRs. In such a scenario, the largest grains generated by planetes- 
imal collisions would remain in resonance, while the smaller grains would feel the effects of 
radiation pressure and immediately occupy orbits with larger semi-major axes and different 
eccentricities. Depending on the collision rate of the small grains, PR drag may act quickly 
enough to bring these small grains back into resonance, possibly creating resonant structures 
with even larger transit signals than we modeled. The lifetime of a resonant planetesimal 
belt would be even shorter than the lifetime of our modeled dynamically cold planetesi- 
mal belt, thou gh, since the mean mot ion resonances would enhance the collision rate of the 
planetesimals (jStark fc Kuchnerll2009l ). eroding the disk at a faster rate. 



4.2. Photometric detection & sources of confusion 

Kepler was designed to detect the transit signature of an Earth-sized planet, requiring 
precise photometry on the order of one part in 10 million, but only over timescales of a 
few hours. By modeling and removing the long-term systematic variations of Kepler, phase- 
folding a light curve obtained over many orbits, and binning the data to increase the signal to 



noise ratio, one can improve the long-term photometric precision. For example, IWelsh et al. 
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( 120 lOl ) detected the exoplanet-induced ellipsoidal stellar variability of Hat-P-7, with a signal 
amplitude ~ 10~ 4 , after 33.5 days of observation, or 15 orbits of the perturbing exoplanet 
Hat-P-7b. 

Our models of exozodiacal resonant ring structures have orbital periods much longer 
than that of Hat-P-7b; a planet orbiting at 0.1 AU from a Sun-like star takes 11.6 days to 
complete an orbit. Since the dust sublimation radius prohibits us from considering planets 
on smaller semi-major axes, the detection of an exozodiacal resonant ring structure around 
a Sun-like star with Kepler requires a deep understanding and precise modeling of Kepler's 
long-term systematic variability. Such modeling is beyond the scope of this paper, as the 
Kepler Pre-search Data Conditioning (PDC) software has not been made public. 

In addition to modeling the long-term stability of Kepler, future efforts to uncover the 
transit signal of an exozodiacal resonant ring within the Kepler data will need to contend 
with any possible sources of confusion. For example, in rare circumstances stellar variability 
from star spots may resemble our exozodi transit signals if the stellar rotation period is close 
to one half or any integer multiple of the planet's orbital period and the distribution of star 
spots creates two minima that lead and trail the planet. Although there is evidence for 
the synchronization between a star's rotational period and a Jupiter-mass planet's orbital 



period near 0.1 AU ( lLanza et al.ll2009f ). differential rotation of the star spots as a function of 
stellar latitude and their relatively short lifetimes make this scenario short-lived and therefore 
distinguishable from the fixed-period exozodi signal. 

Additional planets may also affect observational data to create seemingly broad transit 
curves under special circumstances. For a system with two planets in the 2 : 1 MMR, the 
phase folded light curve of the inner planet would also show the transit of the exterior planet. 
The exterior planet's transit would consistently occur at r oughly the same phas e of the inner 



planet and would not produce a broad transit signal. But iLissauer et al.l (120111 ) showed that 
many multiple planet candidate systems detected by Kepler exhibit planet candidates that 
are near, but not within, mean motion resonant configurations. For such a system, the phase 
folded light curve of an inner planet would also show the transit signatures of all exterior 
planets near j : 1 resonances, but those transits would drift in planetary phase. This drift 
would broaden the transit of the exterior planet by 
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in the limit \jT{ — T Q \ <C T« <C i b s , where t ^ s is the total observation time, Tj is the 
orbital period of the inner planet selected for phase-folding, and T Q is the orbital period 
of any exterior planet near a j : 1 mean motion resonance. If two roughly terrestrial-sized 
exterior planets near j : 1 resonances transit ~ 90° before and after the inner planet's transit, 
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the composite broadened transit signatures may look similar to our resonant ring transit. A 
planet near an interior 2 : 1 mean motion resonance could also contribute a broadened transit 
curve that occurs twice per folded orbit, similar to our resonant ring transit curve. In these 
cases, a simple periodogram and the non-phase folded light curve should distinguish between 
additional planets and a resonant ring structure. 

As Figure [3] shows, the transit curves of resonant rings exhibit several features that may 
aid in their disambiguation, including minima that may be separated by more than 180° in 
planetary phase and a broad maximum near a planetary phase of 180° that peaks at a value 
less than the overall maximum. However, low signal to noise may make identifying these 
features difficult. Future missions with better photometric precision than Kepler over longer 
time scales would facilitate the detection of these disk features. Simple two- or three-band 
transit photometry may also help identify the non-gray transit features expected for small 
dust grains, e.g. the predicted color asymmetry between the leading and trailing clumps as 
discussed in Section! 



Because resonantly trapped dust may originate from parent bodies with a distribution 
of inclinations up to a few tens of degrees, the transit signature of a resonant ring structure 
may be detectable over a slightly wider range of system inclinations than planetary transits. 
Identifying the signal unambiguously without any accompanying planetary data may be a 
difficult task for Kepler. Detection of the weak photometric variations of the disk (~ 10~ 5 ) 
over orbital time scales will require multiple rotations of the ring (i.e. phase-folding) and 
broad time binning to achieve adequate SNR. Prior knowledge of the planets orbital period 
provides us with the proper phase-folding period to perform this task. Additionally, the 
phase between the disk and planetary transits provides the primary evidence linking the 
observed signal to a resonant ring of dust produced by a planet. 

Detection of disk transit signals with Kepler may open up a new region of parameter 
space for exozodi detection. Spitzer was sensitive only to den se exozodis (~ 1400 zodis) 



within 1 AU of nearby solar-type stars (jBeichman et al.l I2Q06I ) . Within the decade, the 



James Webb Space Telescope may extend this detection limit to stars as distant as the typical 
Kepler field star, assuming similar spectral calibration uncertainties. Within the next few 
years, the Large Binocular Telescope I nterferometer may de tect exozodiacal disks as tenuous 



as 10 zodis, but only for nearby stars (ILawson et al.l 120091 ). In the immediate term, Kepler 



observations of disk transits provide the only method for probing asymmetric structure in 
exozodiacal clouds around distant stars, possibly enabling the detection of exozodiacal dust 
disks with densities ~ 100 zodis or less. 
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4.3. Caveats 



The signal amplitudes shown in Figure |6] assume blackbody dust grains (Q a bs = 1) 
with an albedo of cu = 0.5. Dust in the zodiacal cloud has an albedo significantly lower, 



~ 0.1 (IDumont fc Levasseur-Regourdlll988l ; lHahn et al.ll2002l ). while rece nt observations of 
the outer regions of debris dis ks reveal dust albedos ranging from ~ 0.05 (IKrist et al.l 12010 



Kalas et all 120051 1 to > 0.4 ([Weinberger et al.l Il999l : iDebes et al.l 120081 ). As long as our 
assumption of Qabs = 1 is approximately valid at visible wavelengths, an albedo of would 
only reduce the modeled signal by a factor of 2. 

As previously mentioned, our treatment of grain-grain collisions ignores grain fragmen- 
tation. We expect fragments with /3 > 0.5 to affect the results only marginally, since these 
grains leave the system on dynamical timescales much shorter than the collision time. To 
first order, including fragments with < 0.5 would increase the number of small grains 
in the inner disk and decrease the number of large grains, since the additional fragments 
would enhance the collision rate. As a result, massive planets that readily trap small grains 
(M p ~ 1 Mj) may produce an even more pronounced ring structure, while planets that do 
not trap such small grains should exhibit a weaker ring structure. 



The effects of fragments with (3 f» 0.5 is less clear. IStrubbe fc Chiang! (120061 ) showed 
that this population of highly eccentric, barely-bound grains is unusually long-lived compared 
with other grains and dominates a disk's optical depth exterior to its birth ring. Production 
of these smallest fragments in the inner disk likely produces results similar to the /3 < 
0.5 fragments, but their near- unity eccentricities would make resonant trapping less likely. 
Future work to improve grain fragmentation in the collisional grooming algorithm will shed 
light on these issues. 

Perhaps the largest caveat of our models is that we assume a single planet on a circular 
orbit co-planar with the disk. A planet with non-zero eccentricity or inclination will trap dust 
less efficiently than its circular, un- inclined counter p art an d/or form resonant structures with 
less asymmetry. For example, iDeller fc Maddisonl ( 120051 ) showed that Jupiter-mass planets 
orbiting at 45 AU with eccentricity < 0.1 trap dust into clumpy resonant structures similar to 
that shown in the top- left panel of Figure [H but the same planet with an orbital eccentricity 
equal to 0.3 produces a ring without clumps. Massive planets on eccentric orbits may also 
create resonant dust rings that rotate at half the rate of the planet's o rbital frequency, re- 
quirin g two orbits of the planet to observe both exozodi transit minima (IKuchner fc Holman 
2003h . 



Additional planets may also reduce the population of dust in a planet's resonances. 
Gravitational perturbations by nearby planets could reduce a planet's dust trapping effi- 
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ciency or reduce the amount of time dust spends in resonance. Additional planets located 
between the inner planet and the planetesimal belt could also eject in-spiraling dust grains 
before they reach the inner planet's MMRs. Jupiter and Saturn shield the inn er few AU 
of our Solar System, ejecting more than 80% of in-spiraling Kuiper Belt dust (ILiou et al. 
19961 ; iMoro-Martm fc Malhotral 120031 ; iKuchner fc Starkl |2010| ). Planetary system architec- 
tures such as this would greatly reduce the likelihood of a massive exozodiacal cloud being 
fed by a distant Kuiper Belt- like source of dust. 

Our models also assume that 100% of the dust is generated by a dynamically cold pop- 
ulation of planetesimals. Dust generated by a dynamically cold population of planetesimals, 
like the asteroid belt, has a higher chance of being resonantly trapped by a planet when 
compared to dust from a dynamically hot population, like an inclined, eccentric population 
of comets. In the zodiacal cloud, the relative amounts of dust at 1 AU supplied by asteroids 
and comets is unknown; estimates of the relative dust contribution from th e asteroid belt 
range widely from ~ 5% to ~ 50% (jlpatov et al.ll2008l ; iNesvorny et al.l 120101 ) . However, we 
predict that exozodiacal transit signals will only be detectable in the near future in systems 
that are by definition much different from the Solar System, with Jupiter-mass planets or- 
biting interior to the orbit of Mercury. We should not expect the planetesimal distributions 
in these systems to mirror our own. 



5. Conclusions 

We used a collisional grooming algorithm to model the dust distributions of collisional 
exozodiacal clouds perturbed by single planets on circular orbits around Sun-like stars. We 
imaged these dust disks edge-on at visible wavelengths as a function of planetary phase 
to synthesize transit light curves of structured exozodiacal clouds. The "clumps" in an 
exozodiacal cloud, formed by the planet's resonant trapping of in-spiraling dust, occult 
starlight and typically create two broad transit minima that lead and trail the planetary 
transit. 

We measured the transit depth of these light curves as a function of planet mass and 
semi-major axis, disk density, and the dust grains' optical scattering parameters. We found 
that Jupiter-mass planets produce the largest transit signals for exozodi structures, with am- 
plitudes up to ~ 10~ 4 . Neptune-mass planets and super-Earths can create exozodi transit 
amplitudes ~ 10~ 5 and ~ 10~ 6 , respectively, while Earth-mass planets produce weak reso- 
nant ring structures for zodi levels > 10 zodis. Planets more massive than Jupiter produce 
exozodi transit signals with amplitudes less than 10 -4 . 
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Disks with resonant ring structures possibly detectable with Kepler require optical 
depths approaching ~100 zodis. Compact disks this massive could have lifetimes of 100 
Myr or longer, consistent with a non-transient scenario. The dust could also originate from 
more distant sources and migrate into the inner regions of disks. 

For a Sun-like star, exozodi transit signals vary on timescales greater than a few days. 
Detection of an exozodi transit signal from Kepler data will require accurate modeling of 
Kepler's long-term systematic variability as well as any stellar variability. Future multi-band 
transit photometry with improved precision over timescales of many weeks would aid in the 
detection of the transit signatures of resonant exozodi structures. 

Asymmetric structure in exozodiacal cl ouds may be a critical source of astrop hysical 



noise for future exo-Earth imaging missions (IDefrere et al.ll2009t lLawson et al.l 120091 ). This 



new method of resonant ring transit detection provides the only current probe by which we 
can place upper limits on the degree of exozodi structure. This method may also open up a 
new region of parameter space for exozodi detection around distant stars. 

The author thanks Alycia Weinberger, Marc Kuchner, Evgenya Shkolnik, Ruth Murray- 
Clay, John Debes, and Brian Jackson for helpful discussions. This work was supported by 
the Carnegie Institution of Washington. 
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